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ABSTRACT 

In this Letter, we present the first evidence in the BATSE data for a prompt high-energy (25-300 keV) 
afterglow component from a 7-ray burst (GRB), GRB980923. The event consists of rapid variabilty 
lasting ~ 40 s followed by a smooth power law emission tail lasting ~ 400 s. An abrupt change in spectral 
shape is found when the tail becomes noticeable. Our analysis reveals that the spectral evolution in the 
tail of the burst mimics that of a cooling synchrotron spectrum, similar to the spectral evolution of the 
low-energy afterglows for GRBs. This evidence for a separate emission component is consistent with the 
internal-external shock scenario in the relativistic fireball picture. In particular, it illustrates that the 
external shocks can be generated during the 7-ray emission phase, as in the case of GRB990123. 
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1. INTRODUCTION 

The distinction between 7-ray bursts (GRBs) proper 
and their afterglows may be one between emission from 
shocks internal to the relativistic outflow and an external 
shock with a circumsource medium (Rees and Meszaros 
1994; Sari and Piran 1997; Meszaros and Rees 1997). From 
an observational viewpoint, the relation between the burst 
proper and afterglow emissions is not well understood. 
However, the recent identification of a simultaneous op- 
tical counterpart to GRB990123 made with ROTSE (Ak- 
erloff et al. 1999) has shown that in some cases emissions 
from internal and external shocks can overlap, in agree- 
ment with the work of Sari and Piran (1999a). 

One way to assess this relation has been the compari- 
son of the 1-10 keV X-ray emission during the burst with 
that observed in the afterglow. Backward extrapolations 
of the BeppoSAX X-ray afterglow light curves to the decay 
part of the burst indicate that the late X-ray emission in 
the burst and the afterglow are consistent with a single 
power law decay curve indicating that the tail of the burst 
evolves continuously into the afterglow (Costa et al. 1997; 
Piro et al. 1998; Nicastro 1998; In't Zand 1998). These 
results appear to be in contradiction with the initial dis- 
appearance of X-ray emission in GRB780506, which was 
followed ~ 400 s later by their reappearance, supporting 
the case of afterglow emission disconnected from the burst 
(Connors and Heuter 1998). Elsewhere, evidence for ex- 
tended hard 7-ray emission has been reported for ~ 10% of 
the bursts detected by PHEBUS (Tkachenko et al. 1995). 



OSSE has provided upper limits on high-energy (> 300 
keV) post-burst emission from GRB950421 (Matz et al. 
1996). Analysis of a subset of GRBs dubbed "FREDs" 
(Fast Rise Exponential Decay) from early in the BATSE 
mission has shown a softening trend in the decay of the 
burst (Bhat et al. 1994). More recently, Burenin et al. 
(1999) report GRAN AT/SIGMA observations of a soft 7- 
ray tail lasting ~ 1000 s. 

We present BATSE observations of GRB980923, which 
show that after a period of ~ 40 s of strong variability, the 
signal abruptly enters a phase of smooth decay. During 
this decay phase, which lasts ~ 400 s, the properties of 
the 7-ray emission are very similar to those of afterglows 
at lower energies. Our results support the idea that, in 
this case, the afterglow emission started during the burst, 
consistent with the theoretical prediction of early high- 
energy afterglow by Sari and Piran (1999b). The spectral 
evolution of this 7-ray afterglow indicates that it is the 
early, high-energy part of subsequent X-ray and optical 
afterglow. 



2. GRB980923: BATSE OBSERVATIONS 



On 1998 September 23 at 20:10:47.5 UT BATSE was 
triggered on the 64 ms time scale by a long and intense 
GRB (trigger 7113), shown in Figure 1. The burst was not 
in the field of view of the BeppoSAX WFC. RXTE/PCA 
made a single observation in stare mode in response to 
the BATSE Rapid Burst Response (Kippen et al. 1997) 
system, but observed no change in count rate. A full scan 
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of the GRB error circle by RXTE/PCA was not possible 
at the time. 

The event consists of complex variability lasting <~ 40 s, 
followed by a smooth emission tail. The 50-300 keV peak 
flux of the burst on the 64 ms time scale is (1.16 ± 0.02) x 
10 -5 ergs sec -1 cm~ 2 , ranking twelfth in brightness among 
all BATSE GRBs. The fluence of the burst (> 25 keV) is 
(4.84 ± 0.02) x 10 -4 ergs cm~ 2 , ranking third among all 
BATSE GRBs. The fluence in the tail (> 40 s) is - 7% of 
the total burst fluence. 

3. TEMPORAL AND SPECTRAL ANALYSIS 

Typically GRBs are short enough that prior and post 
background intervals can be chosen so that the behav- 
ior of the background during the burst emission interval 
can be estimated through interpolation. In the case of 
GRB980923, however, the duration of the emission tail is 
unknown, therefore routine background subtraction is not 
applicable. We use as background the average of the count 
rates registered when the spacecraft (CGRO) is positioned 
at the point closest in geomagnetic latitude to that at the 
time of the burst on days before and after the burst trigger, 
as described in Connaughton ct al. (1999). Background 
rates were computed using both CONT (16 channel, 2.048 
s time resolution) and DISCLA (4 channel, 1.024 s time 
resolution) data types for comparison. 

The two detectors which have the highest burst count 
rates also have the smallest source angles to Vela X-l. 
The Vela pulsations were modeled using a 7 term har- 
monic Fourier expansion centered on a barycentric fre- 
quency equal to 3.546 mHz, as tabulated from BATSE pul- 
sar monitoring (M. Finger, private communication). The 
model was fit to the background subtracted DISCLA data 
in channel 1 (25-50 keV) during a pre-burst interval rang- 
ing from —405 to —5 s and a post-burst interval from 300 
to 1000 s, having a x 2 /dof of L12 (1067 dof). Fourier 
spectra of the residual rates in the fitting intervals are 
consistent with that expected from Poisson fluctuations. 
We also modeled the pulsations in the CONT data (used 
for spectral modeling) using the same procedure and time 
intervals. Fits were made in each channel for channel 1 
through 10 with an average reduced x 2 equal to 0.95. 

3.1. Temporal Modeling 

We modeled the tail count rates with a power law (PL) 
of the form A(t — to)' 3 - Separate fits were made in five 
energy ranges (25-50, 50-100, 100-300, 25-300, and > 300 
keV) using the DISCLA count rate data from 40 to 298 s 
after the trigger time. For comparison, fits were also made 
using CONT data, binned in energy to match the DISCLA 
channel boundaries. A series of fits were made in each en- 
ergy range for t = [-14.995,32.109] at 1.024 s intervals 
using a \ 2 minimization algorithm with two free param- 
eters, A and /3, while the value of to was held constant 
for each fit. From the fits in the energy ranges below 300 
keV, we find best-fit (i.e., minimum \ 2 ) values for (3 in the 
range —1.81 to —1.64 (with typical uncertainties ±0.03) 
for to in the range 9.6-12.1 s. For the high-energy (> 300 
keV) emission we find (3 — — 1.85 ± 0.53 for to = 10.1 s 
(x 2 /dof = 0.2). The 99% confidence intervals in to range 
from —4 to +18 s, and do not change much with energy. 



The best fit model in the 25-300 keV range using DISCLA 
data gives = -1.81 ± 0.02 (x 2 /dof = 1.22) for t = 9.6 
s (see Figure 2). Both 68.3% (dashed) and 99% (solid) 
confidence intervals for to are shown in Figure 2 with the 
25-300 keV time history and the PL models from the range 
of to values. The PL fits, and the obvious requirement that 
the PL emission cannot exceed the total burst signal, in- 
dicate that the PL decay started ~ 20-25 s after the burst 
trigger. Note that there is no reason why t for the after- 
glow should be 0, i.e. the time of trigger, since the (inter- 
nal shock) GRB emission and (external shock) afterglow 
are physically different and thus somewhat independent in 
timing. 

3.2. Spectral Modeling 

For the spectral modeling we used the Vela-corrected 
CONT data (22-1880 keV) and constructed a single time- 
integrated spectrum for a time interval spanning 40.813 
to 102.253 s in the tail. At the time of observation, the 
onboard LAD look-up table was in soft mode, mapping 
low energies into higher CONT channels and thus provid- 
ing suitable statistics below 50 keV. To examine the spec- 
tral evolution in the tail we generated ten spectra with 
8.192 s time resolution (four 2.048-second CONT bins) 
from 39.789 to 121.709 s in the tail. Although the tail of 
the burst appears to last beyond ^250 s, we selected our 
source intervals such that the model fit parameters are 
reasonably constrained. The fits were made using the for- 
ward folding technique. The photon model that best rep- 
resents the data was determined by comparing fits made 
with a single PL (2 free parameters) and a smoothly bro- 
ken power law (SBPL) (4 free parameters) using the Ax 2 
statistic (Band et al. 1997). We obtain Ax 2 = 67 with a 
chance probability of obtaining a value equal to or greater 
than Ax 2 of - 10~ 15 . Therefore, we used the SBPL to 
model the spectra for GRB980923. 

The parameters of the SBPL are the low- and high- 
energy photon indices ai ow and «high, and the break en- 
ergy Eb. From the time-integrated fit, we obtain ot\ ow = 
-1.66±0.04, a high = -2.20i0.ll, and E h = 203i59 keV 
(x 2 /dof = 1.06). The GRB function (Band et al. 1993) 
was also fit (x 2 /dof = 1.07) giving nearly identical param- 
eter values (a G RB = -1.59 i 0.05, /3 G rb = -2.192 i 0.14, 
E p = 271 i 52), similar to values commonly found in GRB 
spectra (Preece et al. 1999). However, the GRB function 
proved to be less robust than the SBPL in the time re- 
solved fits. From the ten fits made in the tail, we find that 
aiow and ahigh do not vary significantly, with weighted av- 
erage values —1.68 i 0.03 and —2.09 i 0.07, respectively. 
The break energy was held fixed in several of these fits to 
allow meaningful estimates of the photon indices. There- 
fore, we repeated the ten fits keeping ai ow and ahigh con- 
stant at their average values, allowing only the break en- 
ergy and amplitude as free parameters. decreases with 
time (see Figure 4), indicating that the spectrum retains 
a constant shape as it shifts in time to lower energies. 

Time resolved spectra (2.048 s resolution) during the 
main burst emission were also modeled with a SBPL us- 
ing CONT data corrected for deadtime effects and Vela 
pulsations. Figure 3 shows the evolution of ai ow and 
ahigh during the main burst and the tail. An abrupt 
change in the spectral form is quite evident where the 
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main burst emission drops sharply and the tail of the burst 
becomes detectable at ~ 40 s. The time integrated spec- 
trum of the main burst gives GRB function parameters 
( a GKB = -0.61±0.01, /3 G rb = -2.95±0.08, E p = 364±3), 
typical for GRBs (Preece et al. 1999). 

Careful inspection of Figure 3 reveals a softening in the 
low-energy index, ai ow , near 14 s, where the burst inten- 
sity drops to roughly the same level as observed in the 
tail (see Figure 1 and 2). Beyond ~ 14 s «i ow appears to 
maintain slightly steeper values, suggesting an overabun- 
dance of low-energy flux that was not present (or at least 
very weak) at the time of the burst trigger, which may be 
attributed to the rise of the afterglow since the decay must 
begin ~ 20-25 s after the trigger time. 

4. DISCUSSION 

In the external shock scenario of the relativistic fire- 
ball model, electrons arc accelerated to a PL energy dis- 
tribution N(-f c )dj c oc 7c~ p c?7c by a forward shock gener- 
ated when the relativistic shell encounters the surround- 
ing medium, where j c is the electron Lorentz factor. The 
synchrotron spectrum of a relativistic shock with a PL 
electron distribution can be described by four PL regions 
segmented according to v a < v c < v m in the fast-cooling 
regime, and v a < v m < v c in the slow-cooling regime (Sari 
et al. 1998), where v a is the self-absorption frequency, v c 
is the cooling frequency, and ^ m is the characteristic syn- 
chrotron frequency. For adiabatic evolution, v c evolves as 
t- 1 / 2 , while u m oc i~ 3 / 2 (Sari et al. 1996). In the fast- 
cooling regime, the evolution may be radiative, in which 
case oc t~ 12 / 7 . Adopting the notation of Sari et al. 
(1998) and ignoring self-absorption, the observed spectral 
flux in the fast-cooling regime is given by 

p _ f ^,max(zV^c)~ 1/2 , V C <V < U m , 

(1) 

Similarly, the flux in the slow-cooling regime can be writ- 
ten as 

F = { ^,ma X (^/fm) _(P_1)/2 , V m < V < V c , 

" I F Vjmax K/^m)~ (P ~ 1)/2 0V^c)" P/2 , U>U C . 

(2) 

Three characteristics of the synchrotron spectrum allow us 
to distinguish between the fast- and slow-cooling regimes: 

1) the change in spectral slope across the break frequency, 

2) the time dependence of the break frequency, and 3) the 
relation between the temporal and spectral PL indices. 

For slow-cooling we expect a change in the spectral slope 
across u c equal to p/2 — (p — l)/2 = 0.5; for fast-cooling, 
the change across v m equals p/2 — 1/2. The value of p can 
be obtained from p/2 = — (cthigh + 1), leading to p = 2.4 ± 
0.11, typical for afterglows. The expected change in slope 
in the fast-cooling regime then equals 0.7 ± 0.11. From 
the time- integrated spectral fit, we find a slope change 
of 0.54 ± 0.12; from the time-resolved fits we obtain an 
average change in slope of 0.42 ± 0.09. These values are 
in good agreement with slow-cooling, but do not exclude 
fast-cooling very strongly. 



If we are observing a cooling break then we can expect 
the break energy to evolve as t -1 / 2 (Sari et al. 1996) 
for adiabatic evolution. The time-dependent break energy 
was modeled with a PL of the form E (t — t n ) s in the same 
manner as the decay of the intensity. We find a PL index 
S = -0.52 ± 0.12 (x 2 /dof = 1.12) 'for t = 32.109 s (see 
Figure 4) ; slightly beyond the 99% confidence interval, but 
within a 99.9% interval and consistent with the PL decay 
starting 20 s after the trigger time. To obtain a fit with 
S = —3/2 or S = —12/7, corresponding to the two possible 
behaviors of v m > v c (namely the adiabatic and radiative 
cases, respectively), would require to to be negative, and 
is thus excluded by the data. We conclude that the break 
energy corresponds to the cooling frequency, i.e. the tail 
is in the slow-cooling regime. 

If we insert the time dependence of v c and v m into the 
expressions for F v and let a = (p — l)/2 and a 1 = p/2, 
where a and a' are the low- and high-energy PL indices 

of the energy spectrum, then we have F v oc v~ a t 2 Q for 

/ _3 / , 1 

v <v c and F v oc v a t 2 Q +2 for v > v c . Therefore, the 
relationships between the temporal and spectral indices 
are (5 = -3a/ 2 for v < u c and j3 = -3a' / 2 + 1/2 = 
— 3a/2 — 1/4 for v > v c . The values of ai ow and ahigh from 
the time integrated spectral fits give an expected values of 
/? = — 0.99±0.04 for the low frequency light curve (y < v c ), 
and (3 = — 1.3 ±0.11 for the high frequency light curve [y > 
v c ). Similar values are derived when the average spectral 
indices are used. Since the PL decay must begin later 
than 20 s after the burst trigger, we can expect the decay 
index to be shallower than —1.64. For to = 32.109 s we 
obtain /3 = -0.98±0.02 (25-50 keV), in agreement with the 
value computed from the spectral index and in accordance 
with the value of t that gives the correct break energy 
PL index. Although this fit does not have the lowest x 2 
(% 2 /dof = 1.54), it does satisfy the requirement that the 
PL can not exceed the total burst count rate. The same 
scenario holds for the high frequency light curve. Our 
measured value (3 = —1.85 ± 0.53 in the high frequency 
regime (> 300 keV) is consistent with (3 = —1.3 derived 
from the spectral index. 

In the fast-cooling regime, two cases must be con- 
sidered. For adiabatic blast wave evolution, we expect 
F v oc t( 2 ~ 3 ?>)/ 4 = t -1 - 3 above v m which agrees with the 
measured high-energy decay. However, the expected be- 
havior below v m is F v oc t~ 1//4 , which strongly disagrees 
with the measured 25-50 keV PL index by at least 0.76 
for all values of to- For radiative blast waves, the behav- 
ior above v m is F oc t( 2 ~ 6 P)/ 7 = t -1 - 8 , in agreement with 
the data; however, below v m we expect F oc t~ 4 / 7 , which 
also strongly disagrees with the data for all values of to- 
la summary, we conclude that fast-cooling scenarios with 
v c < v m do not agree with the evolution of the tail of this 
burst. 



5. CONCLUSION 

We find similarities between the properties of 
GRB980923 and those of GRB920723 (Burenin et al. 
1999). Both exhibit long soft 7-ray tails, however the 
tail of GRB920723 last longer by a factor 2 and decays 
less steeply (fi — —0.7). Burenin et al. report a change 
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in spectral index equal to ~ 0.8, which is strikingly sim- 
ilar to the abrupt change we observe in a\ OV! . From our 
spectral analysis we have been able to identify which blast 
wave evolution regime is favored in the production of the 
high-energy tail. 

These results suggest a behavior not unlike that of af- 
terglows observed at lower energies. In particular, the evo- 
lution of the spectrum in the tail of the burst appears to 
mimic the evolution of a synchrotron cooling break in the 
slow-cooling regime. This implies that the transition from 
fast to slow-cooling (i.e., when v m = v c ) can take place 



on (short) time scales comparable to the duration of the 
burst. If the synchrotron emission from external shocks is 
the source of the "classical" afterglow, then our analysis 
provides evidence that, at least in some bursts, the after- 
glow may begin during the 7-ray emission phase and may 
reveal itself in the GRB light curves. 
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Fig. 1.— DISCLA time history of GRB980923 (> 25 
keV). 
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Fig. 2.— DISCLA time history of GRB980923 (25-300 
keV) plotted logarithmically with a range of PL mod- 
els extrapolated backwards in time through the burst 
proper. The confidence intervals for to are designated 
by the dashed (68%) and solid (99%) curves. The dot- 
ted line is the fit for t = 9.581 s (j3 = -1.81 ± 0.02, 
X 2 /dof-1.22). 
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Fig. 3. — SBPL low- and high-energy photon indices (with 
one-sigma uncertainties) as a function of time since trig- 
ger. 
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Fig. 4. — Break energy as a function of time since to on 
logarithmic axes. The dashed line is the theoretical slope 
of —1/2 for synchrotron cooling. The solid line is the fitted 
PL slope -0.52 ± 0.12 (t = 32.109, x 2 /dof = 1.12). 



